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The electronic structure of pyramidal shaped InAs/GaAs quantum dots is calculated using an 
eight-band strain dependent k • p Hamiltonian. The influence of strain on band energies and the 
conduction-band effective mass are examined. Single particle bound-state energies and exciton 
binding energies are computed as functions of island size. The eight-band results are compared with 
those for one, four and six bands, and with results from a one-band approximation in which m e ff{f) 
is determined by the local value of the strain. The eight-band model predicts a lower ground state 
energy and a larger number of excited states than the other approximations. 



I. INTRODUCTION 



Semiconductor quantum dots made by Stranski-Krastanow growth have been of great interest over the past few 
years. Such heterostructures are made by epitaxially depositing semiconductor onto a substrate of lattice mismatched 
material. The deposited material spontaneously forms nm-scale islands which are subsequently covered by deposition 
of the substrate material. In this way electrons and holes may be confined within a quantum dot of size 10 nm or loss. 
The islands have a pyramidal shape with simple crystal planes for their surfaces. The presence of strain significantly 
alters the electronic structure of the quantum dot states. Theoretical studies of strained islands have employed various 
degrees of approximation to the geometry, strain distribution, and electron dynamics, ranging from single-band.models 
of hydrostatically strained islands, to multiband models including realistic shapes and strain distributions JjlJ 

In this paper we consider an InAs island surrounded by GaAs. Due to the large lattice mismatch ( w 7%) the 
strain effects are substantial. The influence of strain is compounded by the fact that InAs has a narrow band gap 
{Eg = 0.418 eV), implying strong coupling between the valence and conduction bands. This provides a compelling 
reascm to use an eight-band model. To date, strain-dependent eight-band calculations have been done for quantum 
wiresa, but not for dots. 




FIG. 1. Island geometry. The island geometry is parameterized by the length of the base, b. 

We assume the island is a simple square-based pyramid with 101-type planes for the sides, as shown in Fig. 1. The 
size of the island is parameterized by the length of the base, b. The choice of island shape is somewhat arbitrary. There 
is no clear consensus on the exact shape, and it may vary with the details of growth conditions. Xhe simple pyramidal 
geometry of Fig. 1 was chosen primarily because it has been used in previous calculations,™ hence facilitating 
comparisons. 

An unavoidable consequence of Stranski-Krastanow island formation is that 1 — 2 monolayers of island material 
remains on the substrate surface. This wetting layer is omitted from the calculations because it may be accounted for 
separately. The strain is insensitive to the wetting layer primarily because it is so thin, and also because it is biaxially 
strained to match the substrate lattice. The wetting layer does play a role in the electronic structure since it provides 
a quantum well state that is coupled to the quantum dot state. However, we are most interested in the tightly bound 
quantum dot states. For these states the wetting layer and quantum dot may be treated separately, simplifying the 
analysis of different wetting layer thicknesses. 
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After briefly outlining the calculational methods, we examine the strain-induced band structure, the single particle 
energies as a function of island size, and the exciton binding energy. Finally we compare the eight-band results with 
calculations using one, four, and six bands. 



II. CALCULATION 



The technique used to obtain the electronic structure has been described previously, where is was used for a six-band 
calculation of InP islands embedded in GalnPD Here we will focus on the differences due to the use of eight bands. The 
entire calculation is done on a cubic grid with periodic boundary conditions. First, the strain is calculated using linear 
continuum elastic theory. The strain energy for the systemQ is computed using a finite differencing approximation, 
and then minimized using the conjugate gradient algorithm. 

The electronic [Structure is solved in the envelope approximation using an eight-band strain-dependent k • p Hamil- 
tonian, + H s a The kinetic piece of the Hamiltonian is 
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Po is the coupling between the conduction and valence bands, E c and E v are the (unstrained) conduction and valence- 
band energies respectively, and A is the spin orbit splitting. The 7^'s are modified Luttinger parameters, defined in 
terms of the usual Luttinger parameters, 7^, by 
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where E g = E c — E v is the energy gap, and E p = 27?t,oPq /Ti 

The strain enters through a matrix-valued potential that couples the various components, 
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6jj is the strain tensor, b and (i are the shear deformation potentials, a v is the hydrostatic valence-band deformation 
potential, and a c is the conduction-band deformation potential. 

In addition to the explicit strain dependence in H s , there is a small piezoelectric effect which is included. The strain 
induced polarization_of the material contributes an additional electrostatic potential which breaks the C4 symmetry 
of the islands to C2 .□ 

The energies and wave functions are computed by replacing derivatives with differences on the same cubic grid used 
for the strain calculation. The material parameters and strain in Eq.'s 1-5 vary from site to site. The Hamiltonian is 
then a sparse matrix which is easily diagonalized using the Lanczos algorithm. The calculation is further simplified 
by eliminating unnecessary barrier material, since bound states fall off exponentially within the barrier. 



III. MATERIAL PARAMETERS 



The values used for the various material parameters are given in Table I. All the parameters were set to the 
values corresponding to the local composition, except for the dielectric constant, cr, which was set to the value 
for InAs throughout the structure. Most parameters were taken from direct measurements, however a few merit 
comment. Neither m c nor a v have been directly measured, although the gap deformation potential a g = a v + a c 
has been measured.Q Using the fact that for most III-IV semiconductors a c /a v « 0.1 Jl3 a c and a g can be estimated. 
Another important parameter is the unstrained valence-band offset, E v b defined as -E"„(InAs) — ^(GaAs) in the 
absence of strain. The value used is phased on transition- metal impurity spectra, and is in agreement with the 
value from Au Schottky barrier data.Eil The value is derived using the fact that transition metal impurities are 
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empirically observed to have energy levels fixed with respect to the vacuum, relatively independent of their host 
environment. Thus, by comparing band edges referenced to the impurity levels in two different materials one deduces 
the relative band offsets if the strain could be turned off. The ground state energies of Mn impurities are 0.028 eV and 
0.113 eV above the valence band in InAs and GaAs respectively^, so the InAs valence band is 85 meV above GaAs. 



TART,F, T 

_1_ iT. J_J J_J J_J ±. 


Material parameters. Unless otherwise noted, values are taken from Landolt-Bornstein. 


ei4 is the piezoelectric 




is the relative dielectric constant, the C's are the elastic constants, and a is the lattice constant. 


Parameter 


InAs 


GaAs 


it 


19.67 


6.85 




8.37 


2.1 


ik 


9.29 


2.9 


E,j 


0.418 eV 


1.519 eV 


A 


0.38 eV 


0.33 eV 


E p 


22.2 eV 


25.7 eV 


a g 


-6.0 eV 


-8.6 eV a 


a c 


-6.66 eV 


-9.3 eV a 


a v 


0.66 eV 


0.7 eV a 


b 


-1.8 eV 


-2.0 eV 


d 


-3.6 eV 


-5.4 eV 


(-14 


0.045 C/m 2 b 


0.159 C/m 2 b 


(-B 


15.15 


15.15 c 


Cxxxx 


8.329 x 10 11 dyne/cm 2 


12.11 x 10 11 dyne/cm 2 


Cxxyy 


4.526 x 10 11 dyne/cm 2 


5.48 x 10 11 dyne/cm 2 


Cxyxy 


3.959 x 10 11 dyne/cm 2 


6.04 x 10 11 dyne/cm 2 


a 


0.60583 nm 


0.56532 nm 


E v bo 


85 meV d 





a Reference =1 

b Referenced 2 ! 

c Value for GaAs used. 

d see text. 



IV. BAND STRUCTURE 



Some insight may be obtained by examining the strain induced modification to the band structure. Fig. 2 shows 
the band energies computed using the local value of the strain (i.e. the eigenvalues of Eq. [| with k = 0). Since the 
coupling between conduction and valence bands is proportional to k, for k — the model reduces to a six-band model 
with a decoupled conduction band. The bands are shown for an island with b = 10 nm. The band diagrams for a 
different sized island are obtained by simply rescaling the x-axes. The dominant effect of the strain is that the island 
experiences a large increase in its band-gap due to the considerable hydrostatic pressure. The conduction band still 
has a potential well 0.4 eV deep at the base of the island, tapering to 0.27 eV at the tip. 

The valence band has a more complex structure. If we could somehow turn off the strain, the holes would be 
confined to the InAs by only E v bo = 85 meV. Strain alters this considerably, and makes the dominant contribution 
to the hole confinement. The most notable features are that the valence band is peaked near the tip of the island, 
with another high point near the base and a band crossing in between. This is most clearly seen in the plot of band 
energies along the 001 direction. The plot along the 100 direction near the base (Fig. 2b) shows that there is a slight 
peak in the valence band at the edge of the base, a feature shared with InP/GalnP islands. InP/GalnP islands also 
have such peaks in the valence band ;r and they are sufficiently strong that holes are localized near the peaks rather 
than spread out over the whole island.u In InP/GalnP islands there is also a valence band peak in the barrier material 
above the island which provides a separate pocket in which holes may be confined. From Fig. 2b we see that InAs 
islands have an elevation of the valence band above the island, but inside the island the valence-band edge is even 
higher. Hence we do not expect holes to be trapped in the barrier material. 
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FIG. 2. Band structure based on the local value of the strain, (a) Bands along the 001 direction, through the center of the 
island, (b) Bands Along the 100 direction, through the base of the island. 



unstrained 
m eff =0.023m 




t i . i i i i i i i i i i i i i i i i i . I 

0.03 0.04 0.05 0.06 

Fig. 3 <m e ff>/m 

FIG. 3. Histogram of the conduction-band effective mass within the island. The effective mass is averaged over direction. 



Strain-dependent effective masses may be found by computing the dispersion relation using the local value of 
the strain. Since m e // is anisotropic, it is necessary to average over directions. The valence-band anisotropics are 
sufficiently strong that a hole effective mass is of dubious value. The conduction-band anisotropy is considerably 
smaller, however, making an electron isotropic effective mass a reasonable approximation. Fig. 3 shows a histogram 
of m e ff within the island. Due to the large hydrostatic strain, m e ff is doubled throughout much of the island, 
although there is considerable variation. 



V. BOUND STATES 



The bound state energies were computed as a function of island size using the full eight-band Hamiltonian (Fig. 
4). Because the calculations were performed assuming no wetting layer there are states right up to the GaAs band 
edges. The energies for one and two-monolayer wetting layers are also shown in Fig. 4 for comparison. These were 
calculated as independent quantum wells using the envelope approximation. It should be noted that there will be 
a bound quantum dot state regardless of island size. It is a well known fact that in one and two dimensions an 
arbitrarily weak attractive potential has at least one bound state. In three dimensions there is no such guarantee, and 
hence it is possible to construct a quantum dot which has no bound states. However the wetting layer forms a quasi 
two-dimensional system with the island acting as an attractive potential. Therefore, we expect at least one wetting 
layer state to be bound to the dot, no matter how small it is. 
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For the conduction-band there are only a few bound states in the island, as shown in Fig. 4a. The number of excited 
states in the quantum dot depends on the island's size and the wetting layer thickness. In order to have an excited 
conduction-band state requires b > 10 nm and b > 12 nm for one and two monolayer wetting layers respectively. The 
first excited state is accompanied by a nearly degenerate state. The splitting between these two states varies from 
2 meV to 6 meV for 10 nm < b < 18 nm. The near degeneracy reflects the C4 symmetry of the square island, with the 
splitting due to the piezoelectric effect. A third excited state appears for b > 13 nm and b > 14 nm for one and two 
monolayers respectively. These limits on b are all upper bounds since the actual quantum dot energies will be lowered 
by the coupling to the wetting layer. In addition to the change in energy with size, the spacings vary as well. The gap 
between the ground state and first excited state varies from 60 meV to 95 meV over the range 10 nm < b < 18 nm. 

The valence-band states are more strongly confined, due to their larger effective mass. Only the first four states 
are shown in Fig. 4b, all of which lie well above the wetting layer energy. The energy spacings vary from a few meV 
to 30 meV over the range of island sizes considered. 

When the island is occupied by an electron and a hole, there will be additional binding energy from the coulomb 
interaction. Ground state exciton energies were computed in the Hartree approximation using eight-band solutions for 
both electrons and holes. That is, the exciton wave function was assumed to be of the form ^ij(r e ,rh) — V*! {r e )ij}^ (fh) ■ 
ip e and ip h were found by self-consistent iteration, with convergence to within 0.1 meV usually taking only two 
iterations. Fig. 5 shows the^exciton binding energy as a function of island size. The results are in good agreement 
with single-band calculationsa, which is not surprising since the coulomb energy depends only on the charge density for 
the electron and hole parts of the wave function. The single-particle electronic Hamiltonian only affects the coulomb 
energy insofar as it alters the charge distribution. The exciton binding energy increases with decreasing island size, 
reaching 27 meV for 6 = 9 nm. Fig. 6 shows the exciton wave function for b = 14 nm. In spite of the complex band 
structure seen in Fig. 2, the electron and hole wave functions appear ordinary. The wave functions are spread out 
over most of the island, with no signs of localization around smaller regions. 




b (nm) 
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FIG. 5. Exciton binding energy versus island size. Exciton binding energies were computed in the Hartree approximation. 




FIG. 6. Electron and hole wave functions for the ground state exciton in the Hartree approximation, with b = 14 nm. 
Surfaces are y^ 8 l^<(r)| 2 equal to 0.1 of the peak value. 



VI. COMPARISON OF APPROXIMATIONS 



There remains the question of whether or not an eight-band model is worth the trouble. Previous authors have 
used various approximations to reduce the number of bands. Grundmann et. al. a treated the electrons and heavy 
holes as single particles moving in the strain induced potential corresponding to the band edges. The effective masses 
were different in the island and barrier materials, but took constant unstrained values within each region. As was 
pointed out by Cusack et. alB the narrow InAs band gap leads to significant band mixing, resulting in large strain 
induced shifts in the effective mass. Based on a pseudopotential calculation, Cusack et. al. set m e // = 0.04mo for 
the conduction band, which is the value predicted for bulk InAs under the average hydrostatic strain in the island. 
The valence-band states were calculated using a four-band model. Note that m e ff = 0.04 m is in good agreement 
with the peak of the distribution for m e ff shown in Fig. 3. 

Unfortunately, a comparison with previous calculations is complicated by the fact that the methods have differed in 
more thafi the Hamiltonian. Different material parameters were used, the strain was calculated differently (continuum 
elasticitya versus valence force field methods), and different numerical techniques were used to solve Schrodinger's 
equation (real space differencingp versus a plane wave basiscl). To more directly compare these different approxima- 
tions, energies were calculated using several different Hamiltonians, but all using the same grid and strain profile. 
For the conduction band the methods considered are (i) setting m e // set to its unstrained values of 0.023mo in the 
InAs, and 0.0665too in the GaAs, (ii) using the value corresponding to the average hydrostatic strain in the InAs 
m e ft = 0.04moi and using m e tf = 0.0665mn in the GaAs, (Hi) using a spatially-varying strain-dependent m e ff[f), 
(iv) using full eight-band Hamiltonian. 

A comparison of the different conduction-band energies for b = 14 nm is shown in Fig. 7a. The dominant feature is 
that the energies decreases as the Hamiltonian includes more physics. For the simple unstrained effective mass only 
a single state is found. With m e // = 0.04too the binding energy of the ground state increases by 30 meV. Also, 
the nearly degenerate first and second excited states are brought below the energy of a one-monolayer wetting layer. 
Using a one-band model with m e //(r) gives energies very close to those for m e ff = 0.04mo. The eight-band results 
are significantly different. Not only is the ground state lower by another 27 meV, but the two nearly degenerate 
excited states are clearly confined. In addition, a third excited state falls below the the one-monolayer energy. The 
one-band models with m e // = 0.04mn and m e ff — m e //(r) both predict E\ — Eq ps 110 meV. The eight-band model 
predicts Ei — E 80 meV. 
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FIG. 7. Confined state energies for an island with b — 14 nm (a) Conduction band computed using constant masses, spatially 
varying mass, and eight-band Hamiltonian. (b) Valence-band states computed using four-, six-, and eight-band Hamiltonians. 

The valence-band energies were calculated using four, six and eight-band models (Fig. 7b). The differences are 
less dramatic than for the conduction-band states. The four and eight-band ground state energies agree to within 
3 meV. For six bands, however, the ground state energy differs from the other two by ~ 40 meV. At first sight this 
is surprising, since one generally expects the more complicated model to produce more accurate results. However, for 
InAs A Eg so if either the conduction or split-off band is to be included, then both should be. The six-band model 
violates this requirement by allowing mixing of the valence-band states, but leaving the conduction band decoupled. 
The eight and four-band models do predict slightly different level spacings, but the basic pattern is the same. E\ — Eq 
is large ( 20 meV for four bands, 16 meV for eight bands) with a smaller spacing between excited states (2 meV and 
8 meV for four and eight bands respectively). 

An interesting way of viewing the model dependence of the energy is to compare with the size dependence. As a 
simple example, comparison of Fig.'s 4a and 7a shows that a b — 14 nm island calculated with m(r) gives the same 
ground state energy as the eight-band model at b — 12.5 nm. Hence, if the uncertainty in the island geometry is 
greater than 1.5 nm we would expect these inaccuracies to dominate the errors in the electronic structure results. 
This gives an indication of the importance of specifying the correct island geometry. 



VII. CONCLUSIONS 



Coupling the valence and conduction bands has a strong impact on the spectrum of InAs quantum dot states. The 
eight-band model gives results significantly different than one, four, and six-band approximations. It predicts larger 
binding energies and strongly confined excited states which do not appear in one-band approximations. 

The results presented here clearly .demonstrate the need for eight-band calculations, or perhaps even more complex 
techniques such as pseudopotentials.Ej While large scale calculations using complex Hamiltonians have become feasible, 
the results are no better than the input parameters used. Accurate agreement between theory and experiment will 
require precise measurements of the island geometry. 

I wish to thank Mats-Erik Pistol, Mark Miller, and Jonas Ohlsson for enlightening comments and discussions. 
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